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^^ ' Abstract 

o 

^ ^ We calculated the spin dependent structure functions g^"" in the valon model rep- 

resentation of liadrons. A simple approach is given for the determination of the twist-3 
^\ , part of 32(2:7 Q^) structure function. Thus, enabling to obtain the full transverse struc- 

. . .' ture function, 52(2;, Q^) of proton, neutron and the deuteron. 



1 INTRODUCTION 

The nucleoli polarized structure functions gi,2{x,Q'^) are important tools in the investi- 
gation of the nucleon substructure. In particular, they are vital elements for the under- 
standing of the spin dependent parton distributions and their correlations. The g2{x, Q'^) 
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structure function is important because it probes both transverse and longitudinally po- 
larized parton distributions inside the nucleon. However, this function is also sensitive to 
the higher twist effects, such as quark gluon correlations which is not easily interpreted in 
pQCD, where such effects are not included and does not fade away even in big Q^ . So it 
seems to be very interesting to calculate it, because it is the only function which is related 
to quark-gluon interaction. Properties of 52 can be clearly established in DIS using the 
Operator Product Expansion(OPE). Therefore, it makes possible to study contributions 
to the nucleon spin structure beyond the simple quark parton model. 
The main purpose of this paper is to calculate transverse spin structure function, (72(2^5 Q^)- 
It requires to consider the twist-2 and the twist-3 contributions. While the twist-2 part is 
well understood, we present a simple method to extract twist-3 part. Since twist-2 part of 
52 requires the knowledge about g\ structure function, a brief review of the latter will also 
be given. In doing so, we utilize the so-called valon model representation of the hadrons. 
Finally we compare our results with the experimental data from E143 [1], E155 [2, 3], and 
also with newly released data from HERMES [4, 5] groups. A comparison of our results 
with other phenomenological models is also given. 

The organization of the paper is as follows. In section 2, we briefly review how to 
calculate the polarized nucleon structure function in the valon model. In section 3 we 
calculate the (72(2;, Q^) spin structure function. We also discuss the numerical results, sum 
rules and the effect of higher twist on 52(2^) Q^)- Then, we will finish with the conclusions. 



2 A brief review of spin structure functions in the valon 
model 

The valon model is a phenomenological model, originally proposed by R. C. Hwa, [6] in 
early 80's. It was improved later by Hwa [7] and Others [8, 9, 10, 11, 12, 13], and extended 



to the polarized cases [14, 15, 16, 17]. In this model a hadron is viewed as three (two) 
constituent quark like objects, called valon. Each valon is defined to be a dressed valence 
quark with its own cloud of sea quarks and gluons. The dressing processes are described 
by QCD. The structure of a valon is resolved at high Q^ . At low Q^, the internal structure 
of the valons can not be resolved and they behave as constituent quarks of the hadron. In 
this model the polarized parton distribution in a polarized hadron is given by: 

5c?f (x, Q") = E /' ^<5GL„ Jy)5gr'^"(^, Q2) (1) 

Jx y y 

where ^G^ualoniv) ^^ ^^^ helicity distribution of the valon in the hosting hadron; in other 
words, it is the probability of finding a polarized valon inside the polarized hadron. The 
Q^ dependence of SG'^^ig^{y) in the next to leading order is very marginal. Moreover, 
it is independent of the nature of the prob. The term (5g™'°"'(x/y, Q^) in Eq. (1) is the 
polarized parton distribution inside a valon. These polarized parton distribution functions 
(PPDFs) are governed by the DGLAP evolution equations. The polarized hadron structure 
functions are obtained via a convolution integral as follows: 

g1{x,Q') = E /' ^<5GLo„(y)ffr"(-,Q') (2) 

valon -^^ y y 
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where (^|"^'°"(^^ q2^ jg ^j^^ polarized structure function of the valon. The details of actual 



calculations are given in [14, 16]. The results obtained for gi{x,Q ) from this model are 
in excellent agreement with all available experimental data. In Fig. 1 we only present a 
sample of the results, and show that how it compares with the more elaborated global fits 
and the other models, we will use those findings here in calculating g2{x, Q^)- 

3 Transverse spin-dependent structure function g2{x,Q'^) 

Polarized Deep Inelastic Scattering(DIS) mediated by photon exchange probes two spin 
structure functions: gi{x,Q'^) and g2{x,Q'^). If the target is transversely polarized, the 
total cross section is a combination of these two structure functions. Transverse spin 
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Figure 1: Polarized proton(deuteron) structure function at Q^ — 5(3)GeV^^ . The results from model 
[14] are compared with the experimental data [18, 19, 20, 21, 22, 23]. 

structure function, g2{x,Q'^), is made up of two components: a twist-2 part, gf^, and a 
mixed twist part, g2{x,Q'^)- It can be written as [24] 



g2{x,Q^)=grix,Q^)+92{x,Q' 



(3) 



where 



Jx oy M y 

The twist-2 part, gf^, comes from OPE. The g2{x,Q'^) receives a contribution from the 



transversely polarized quark distributions hxix^Q"^) as well as a twist-3 component. The 
twist-3 part, Civ^Q^) iii EQ- (4) comes from quark-gluon interactions. These higher twist 
corrections arise from the non-perturbative multi parton interactions, whose contributions 
at low energy increase as tj^, reflecting the confinement. The twist-3 term represents qgq 
correlations. Therefore, any non-zero result for this term at a given Q^ will reflect a 
departure from the non- interacting partonic regime [25]. 
gww jg j-giated to the gi structure function by the Wandzura-Wilczek relation [26] as follows, 



gl^-^{x,Q') = -giix,Q'')+ f\,ix,Q^)^. 

Jo y 



(5) 



There are two important and well known sum rules regarding gi{x,Q'^) and g2{x,Q'^). 
The first one is OPE sum rule: 

r-= f\^gi{x,Q^)dx = ^,n = 0,2,4,... (6) 

Jo ^ 

Yl= f\-g^{x,Q^)dx = l^idn-an),n = 2,4,... (7) 

JO Zn + i 

The second one is Burkhardt-Cottingham sum rule [27], 

/ g2ix,Q^)dx = (8) 

Jo 

The Burkhardt-Cottingham sum rule for g2{x, Q^) at large Q^ was derived from virtual 
Compton scattering dispersion relations. This sum rule does not follow from the OPE since 
the n = sum rule is not defined for g2{x, Q^) in Eq. (7). Since g2 = 32"" ~92, combining 
this with Eq. (5) results in a third sum rule (though not independent). So, we have in 
our disposal the following 3 sum rules 

"\ri^,Q^)dx = (9) 







f\^gr{x,Q^)dx = -\a2 (10) 

Jo o 

I x^g2ix,Q^)dx = ld2 (11) 

Jo o 

where, 02 and ^2 are the matrix elements of twist-2 and twist-3 operators. 

Now we are set to calculate g^" '^{x,Q'^). It will be done in two parts: First we 
evaluate the twist-2 part, g2^{x,Q'^), and then we take up the twist-3 part, g2{x,Q'^). 
Combining the results from the two steps results in the full g2ix, Q'^). 

3.1 Calculation of the twist-2 term, gf'"{x, Q^) 

We begin with Eq. (5). Since gi{x,Q'^) is known in the valon model [9], we utilize those 
results without any additional free parameter and evaluate the twist-2 part of g2{x,Q'^); 
namely g2^{x,Q'^), according to the Eq. (5). The results are shown in Fig. 2 for proton 



and the deuteron. We have included m Fig. 2 the results from other phenomenological 
models for the sake of comparison [28, 29]. The experimental data are from E143 and 
E155 Collaborations [1, 2, 3]. 
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Figure 2: Transverse polarized proton and deuteron structure functions, g2^{x, Q^) at Q^ — 2GeV^ and 
Q2 = 5GeV^. 



3.2 Calculating the twist-3 term, g2{x,Q'^) 

As mentioned before, the function g2{x,Q'^) has two terms: 

-g2{x, Q') = -[' ^{^hT{y, Q') + e(y, Q'))-. 
Jx oy M y 



(12) 



The first term is related to the transverse polarization of quarks in the nucleon and is a 
twist-2 contribution. This term is suppressed by the quark to nucleon mass ratio, so we 
ignore this negligible part and focus on the second term which is the twist-3 contribution 
and related to quark-gluon correlations. 
The moments of g2{x,Q'^) obey the following simple equation [30]: 



^2(n,Q') = L2*o^2(n,Q2) 



(13) 



Where, 



g2{n,Q^)= f\''-'g2ix,Q^)dx, (14) 

2\ 



OsiQj 
L= ,^2\ ^ (15) 

bo = yiVc - \Nf (16) 

7^ = 2iV,(5„_i - i + ^) (17) 

•^n-l = E 7 (18) 

Our task is to solve Eq. (13) with some appropriate initial conditions in the moment 
space. Then we can make a transformation to the momentum space and evaluate the 
twist-3 contribution to the transverse spin structure function. This is done in two steps, 
as is the case in the valon model. The first step involves finding a solution to Eq. (13) 
in a valon. The second step is to convolute the results obtained in the first step with the 
valon distribution in the nucleon. This will give the nucleon structure function. 
Our choice for the initial input function for §2 "^[z, Qq) is : 

-gr'"''iz,Ql) = A6iz-l) (19) 

The justification for this choice is as follows: In the momentum space one can write 

-gr''"'{z,Q') = fm-92(.z,Ql) = f{Q')A5{z - 1) (20) 

Note that for Q^ = Qq we get /(Q^) — ^ f{Qo) = 1 which is apparent from the definition 
of L in Eq. (15). Thus, we arrive at Eq. (19). This simple choice of the initial input for 
gvaion^^^QT-^ stems from the knowledge that it is related to the quark gluon correlations, 
which in turn, is related to the Green function in the momentum space. In the momentum 
space the correlation function is composed of a Dirac Delta term and a function that 
is related to the momentum. Therefore, at the initial Qq we can simply assume that 
gvalon^^^QT-^ is proportional to Dirac Delta function which emphasizes the conservation of 
energy- momentum and the fact that at that low Qq a valon behaves as an object without 



any internal structure. The last point is built in the definition of a valon. So, in the 
moment space the Delta function becomes unity and we can write 



92{n,Ql) = A^ 1 



f211 



all the QCD effects are summarized in A. This coefficient will be extracted from the ex- 
perimental data. With the above choice for the initial input, the moments of §2^ {z, Q^) 
are obtained in the valon, using Eq. (12). An inverse Mellin transformation then takes us 
to the momentum space, giving §2^ {z^Q"^) structure function. This completes the first 
step described above. Meanwhile, a fit to E143 data [1] gives the following functional form 
for 52^'°" (-2, Q^) in the valon, depicted in Fig. 3. 
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Figure 3: The twist-3 part of proton transverse spin structure function in the valon in z space, g2{z, Q^ 



at Q^ = 5Gel/^ 



In Fig. 4 we present x^g^^ for proton and deuteron along with the E143 data [1] data. 



We have also compared our results with those from Bag model [28] and with the results 
obtained by Wakamatsu [31]. 
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Figure 4: The twist-3 part of proton and deuteron transverse spin structure function, x^g2{x,Q^) at 
Q2 = 5GeV^ . 



The second step involves the convolution process which takes us to the hadronic level. 
This is similar to the way we have used to extract g\'^{x,Q'^), and described in a number 
of papers, e.g. see [14, 15, 16]. That is, briefly, to convolute §2^ {z,Q'^) with the valon 
helicity distribution 6G^ '{y) in the corresponding hadron to obtain g2{x,Q'^). After 
adding g2{x,Q^) to the g2^{x,Q'^), the full §2 is obtained and are shown in Fig. 5 for 
proton and deuteron and in Fig. 6 for helium-3 (For calculating the g2^{x,Q'^) , we used 
the polarized proton and neutron helicity distributions which are derived in [32]). 

Finally, in the table 1 we give the numerical values for the twist-2 and twist-3 matrix 
elements. They are defined in Eqs. (9) and (10), respectively. We also include results from 
other sources for the purpose of comparison. The results for the Burkhardt-Cottingham 
sum rule which is defined in Eq. (8) are given in table 2. This integral was evaluated 
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Figure 5: Full transverse polarized proton and deuteron structure function, g2{x,Q^) at Q^ = 5GeV^ 



^He (x,5 GeV^) 



6 - 



4 - 



E142 
Jlab 
Model 



- 



-2 - 



v\i^^^- 



TO" 



10" 



10' 



10° 



Figure 6: Full transverse polarized He'^ structure function, g2{x,Q^) at Q^ = 5GeV^ . We compared our 
results with experimental data from [33, 34] . 
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af 


dP 


ai 


di 


Valon model 


0.0224 


0.0042 


0.010 


0.0037 


MIT bag model [28, 35] 


- 


0.01 


- 


0.005 


QCD sum rule [36] 


- 


-(0.6 ±0.3)10-2 


- 


-0.017 


QCD sum rule [37] 


- 


-(0.3 ±0.3)10-2 


- 


-0.013 


Lattice QCD [38] 


(3 ±0.64)10-2 


-(4.8 ±0.5)10-2 


(13.8 ±5.2)10-'^ 


-0.022 


CM bag model by Song [28] 


0.0210 


0.0174 


0.0087 


0.0067 


E143 [1] 


(2.42 ±0.20)10^2 


(0.54 ±0.5)10-2 


(8.0±0.16)10-'' 


(3.9 ±9.2)10-3 



Table 1: The twist-2 and twist-3 matrix elemet operators ,02 and d2, for the proton and 
the deuteron, calculated in the valon model. Also included the experimental data and the 
results from other theoretical investigations. 





bag model by Song[28] 


E143[l] 


E155[3] 


HERMES 2012[5] 


Valon model 


Jg^^{x,Q^)dx 


-0.0016 


-0.014 ±0.028 


-0.022 ± 0.071 


0.006 ±0.024 ±0.017 


-0.004 


Jgi{x,Q^)dx 


-0.00287 


-0.034 ± 0.082 


0.023 ± 0.044 


- 


0.010 



Table 2: The results for the Burkhardt-Cottingham sum rule. 

in the region of 0.023 < x < 0.9 at Q^ = 5GeV^. We compared our results with those 
experimental data from HERMES in the same region and also with the results from E143 
and E155 in the region of 0.02 < x < 0.8. We compared our results with those from bag 
model too. 

4 Conclusion 

We have used the valon model and calculated the transverse spin structure function for 
nucleon and deuteron. We offer a simple approach for calculating the twist-3 part of 
the transverse spin structure function. It is evident that our results for both the twist-2 
part and for the full transverse spin structure functions are in good agreements with the 
experimental data. We have also provided a comparison of our results with other works. 
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